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ABSTRACT. This paper is an investigation of the stability of some ideal stars. It is in- 
tended as a study in General Relativity, with emphasis on the coupling to matter, eventu- 
ally aimed at a better understanding of very strong gravitational fields and "Black Holes" . 
This contrasts with the usual attitude in astrophysics, where Einstein's equations are in- 
voked as a refinement of classical thermodynamics and newtonian gravity. Our work is 
based on an action principle for the complete system of metric and matter fields, cou- 
pling the metric field to well defined relativistic field models that we hope may represent 
a plausible type of matter. The thermodynamic content must be extracted from the the- 
ory itself. We start with the simplest model of matter and plan to add complications as 
our experience grows. When the flow of matter is irrotational, and described by a scalar 
density, we are led to differential equations that differ little from those of Tolman and his 
followers, but equations derived from an action principle admit a conserved current, as 
well as stronger boundary conditions that affect the matching of the interior solution to 
an external Schwarzschild metric and imply a relation between mass and radius. 

We have been led to propose a complete revision of the treatment of boundary condi- 
tions. An ideal star in our terminology has spherical symmetry and an isentropic equation 
of state, p = ap 1 ', a and 7 piecewise constant. In our first work on this subject it was 
assumed that the density vanishes beyond a finite distance from the origin and that the 
metric is to be matched continuously at the boundary to an exterior Schwarzschild metric. 
But it is difficult to decide what the proper boundary conditions should be and we are 
consequently skeptical of the concept of a fixed boundary. In this paper we investigate 
the double polytrope, characterized by a polytropic index n < 3 in the bulk of the star 
and a value larger than 5 in an outer atmosphere that extends to infinity. It has no fixed 
boundary but a region of critical density where the polytropic index changes from a value 
that is appropriate for the bulk of the star to a value that provides a crude model for 
the atmosphere, with a fairly abrupt transition at the point where the density reaches a 
critical value. All the fields are continuous. The boundary conditions are now natural and 
unambiguous. The existence of a relation between mass and radius is confirmed, as well as 
an upper limit on the mass. The principal conclusion is that all the static configurations 
are stable. There is a solution that fits the sun. The masses of white dwarfs respect the 
Chandrasekhar limit. The application to neutron stars has suprising aspects. 
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1. Introduction 

The problem on which we hope to throw some light in this paper is the application of 
General Relativity to star-like systems that can be described by mass, radius, density and 
pressure; the flow is assumed to be irrotational and to be controlled by a velocity potential. 

In contrast with the traditional treatment we introduce the matter component into 
Einstein's equation for the metric by adding an appropriate matter contribution to the 
Einstein-Hilbert action. The difference, at first sight, seems minor, for the equations 
associated with either approach are nearly identical, especially in the static case. An im- 
portant consequence for the study of equilibrium configurations is that the action principle 
fixes an integration constant that is left free in the phenomenological approach. This re- 
sults in a strengthening of the conditions for matching the interior metric to an exterior 
Schwarzschild metric at the boundary of the star - see Eq.(1.4), or asymptotically at in- 
finity. Consequently, the action principle has some additional predictive power, likely to 
bring it down, perhaps, but worth investigating. 

We do not take it for granted that the matter distribution conforms to the precepts 
of classical thermodynamics; instead it is expected that the interpretation is implied by 
the theory itself. It turns out that a simple choice of interaction does indeed lead to an 
equation of state of the familiar type. The difference in attitude therefore does not, by 
itself, lead to any dramatic contrast with the semi-phenomenological approach. The future 
inclusion of radiative effects may change that. 

The simplest choice of action leads to an equation of state of the form 

p=V, 7=1 + -, (1.1) 
n n 

where a, 7 are constants, eventually piecewise constant. 

The main result of an earlier investigation [F2] was a relation between the mass and 
the radius of any equilibrium configuration. The original purpose of the present paper was 
to study the stability of those configurations. In the course of this work we have become 
somewhat skeptical about the appropriateness of naive boundary conditions; that is, the 
assumed continuity of the metric, regardless of the behaviour of density and pressure that 
it implies, at a fixed boundary. A large part of this paper is directed to a re-evaluation of 
these questions. 

The double polytrope 

Consider the process that leads to the formation of a star, assuming that the primordial 
matter is homogeneous, as is reasonable if stars are a late result of a long process of 
condensation of a hydrogen cloud. Condensation is a result of gravitational attraction 
and the first effect produced by the attraction is an increase in density. All subsequent 
development is ultimately attributable to this primordial increase in the density. If, as 
is always taken for granted in studies of stellar structure, the equation of state is nearly 
isentropic, then the basic, underlying reason for a change in the index n must be the 
variation of density. A popular model is a polytrope with n ~ 3 in a region of moderately 
high density and n > 5 in the outer atmosphere. In the early stages with low density the 



2 



index may be larger than 5 almost everywhere, characteristic of a distribution that extends 
to infinity, but as this would imply a singularity at the origin a change must take place 
near the center. Whatever happens, the primal cause is the variation of density. That is, 
the position of the boundary must be determined by the density, rather than the other 
way around. It follows that, if the index changes abruptly, then it is the result of a rapidly 
changing density, as in 

r i n x {pjp cr ) K + n 2 
" M = ! + („/„„)* ' (L2) 

where K is a suitable large number and p cr is a critical density. (Another representation for 
the approximately piecewise function will be explored at the end, with interesting results.) 

Indeed, in an approximation where the only variables to be taken into account, besides 
the components of the metric field, are density and pressure, this would appear to be the 
only possible approach: the boundary is defined to be the region of critical density. Note 
that this "boundary" need not coincide with the visual boundary of the star. 

In this paper, after attempting a more traditional approach to localizing the surface 
of the star, and remaining unconvinced of the aptness of it, we shall concentrate on trying 
to understand the double polytrope with this type of equation of state. 

The main conclusion is that all the static solutions, with natural boundary conditions 
applied at the center and at infinity, are stable. The white dwarfs respect the Chan- 
drasekhar limit on the mass, not because heavier stars are unstable, but because they do 
not exist. A model for the sun is included without any surprises; the application to neutron 
stars offers new dimensions. 



Summary 

Section 2. An unfamiliar aspect of this work is the use of an action principle for the 
complete system of metric and matter fields. Matter is assumed to be irrotational and 
isentropic, thus fully described by fields of density and pressure, with the action 

A mat ter = J d 4 xV=9~(^ (<7^V, M V> - c 2 ) - V[p\) =: J d 4 x^g C . (1.3) 

The eventual presence of (electromagnetic) radiation will require additional terms to be 
added. The associated energy momentum tensor provides the right hand side of Einstein's 
equation R^ v — g^ v R = 9>txGT^ v . We look for solutions that are spherically symmetric, 
such that the metric, in terms of coordinates t, r, 9 and 0, takes the form 

(ds) 2 = e u {dt) 2 - e x (dr) 2 - r 2 dQ 2 , 

with v and A depending on r and t only. 

Homogeneous poly tropes are characterized by a potential (the internal energy) V[p] 
of the form V[p] = ap 1 , with 7 and a constant, which leads to the equation of state 
p = a(7 — l)p 7 . When 7 is not constant the polytropic equation of state is slightly 
modified in the region of critical density. As is usual, we study the time development of 
the system under the assumption that it is initially in an equilibrium configuration. 
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Our first calculations [F2] postulated a fixed boundary beyond which the density is 
zero and the metric is that of Schwarzschild's exterior solution. The static configurations 
of this model are essentially the same as in the traditional, phenomenological approach, 
except for an important difference with respect to the boundary conditions. In contrast 
with the situation in the usual approach, we must match both of the metric functions v, A 
of the interior solutions to an external Schwarzschild metric. Integration proceeds from 
the center and the boundary is at a point r = R where 



v(R)+\(R)=0. (1.4) 

Within the traditional approach this condition is ineffective since the boundary value of v 
is just an integration constant. 

We wish to calculate the time development to first order in the deviations from equilib- 
rium. This was first done by Chandrasekhar in 1931 [CI], and as far as we know the same 
method has been followed by all later investigators. All these studies are characterized by 
what we think are insufficiently motivated boundary conditions. In the first place it is not 
sure that one knows which of the fields, metric components, density, pressure, should be 
required to be continuous at the boundary. The boundary is not at a fixed point but varies 
from one static configuration to another and with time. Consequently it is unnatural to 
restrict the fluctuations by the condition that the radius remain fixed. Another question 
that imposes itself is that of the mass. We define the mass in terms of the asymptotic 
metric; does it echo the oscillations or does it remain fixed? 

In view of the fact that we have come to view the Schwarzschild solution as the metric 
of a singular, limiting mass distribution [Fl] [F2] [MM], we felt that one way to clarify these 
questions would be to replace the outer Schwarzschild metric with another polytrope, 
with index n > 5 as is appropriate for a mass distribution that extends to infinity. This 
exterior metric rapidly approaches the Schwarzschild metric at moderate distances. But 
the boundary conditions continue to present the same, difficult problems. 

Section 3. We have come to believe that the very idea of a fixed boundary is unnatural, 
and an obstruction to understanding what is going on. For this reason we make a new 
start, with another version of the double polytrope, an ideal star in which the equation 
of state is essentially poly tropic near the center, with an index that is nearly constant, 
but changes to a larger value in a "boundary" region of critical density, and essentially 
constant outside this region. The "surface" of this star is a region in which the index 
makes a sudden or gradual change from one value to the other, as in 

P = a -r+\ P=r, „ = t-zJMf^. (i.5) 

This defines a 2-parameter family of equations of state. For reasons that will be explained, 
the most plausible models have f cr = p cr = 1. 
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Results 

Our earlier calculations were made with a fixed isentropic index for the interior, and 
a matching Schwarzschild metric for the empty exterior. In this case the parameter a can 
be varied by simple rescaling. We investigated these solutions but, having great difficulties 
in selecting the proper boundary condiditions, we abandoned that approach. 

Using the new equation of state, in (1.5), and the natural boundary conditions (regu- 
larity at the center and fall-off at infinity) we re-calculated the static configurations. With 
n = 3 (inside) and n = 6 (outside), solutions were found for values of the parameter a 
ranging from 10 -6 to 1/5.765 and no solutions were found for larger values of a. A relation 
between radius and mass emerges in all cases considered, including: 

(1) Polytropic index 1 < n < 3 (near the center) and 6 (at large distances). 

(2) Index n = 3 and 15. 

All these static solutions appear to be stable to radial perturbations. It may be objected 
that the equation of state used here is somewhat special, not to say ad hoc. We are 
nevertheless justified in concluding that instabilities of polytropes found previously are 
characteristic of a restricted class of boundary conditions; they are not generic. 

An important consequence of the fact that the dynamics is formulated as an action 
principle is the existence of a conserved current. With the usual boundary conditions at 
the center, and natural boundary conditions at infinity, we find that the asymptotic mass 
is a constant of the motion. 

Section 4. An application to the sun predicts the central density and pressure, close 
to the values obtained within the traditional approach. In applications to white dwarfs the 
constant a, a free parameter in other cases, is known. In the case of complete degeneracy 
the polytropic index is equal to 3, and in this case a unique mass is predicted, very close to 
the mass of the sun. This limit, here derived from a theory in which all the stars are stable, 
is close to the limiting value obtained by Chandrasekhar from stability considerations [CI]. 

The instabilities of white dwarfs discovered by Chandrasekhar are difficult to interpret, 
as witness the reservations expressed by Eddington. The first question that comes to mind 
is the future development of a star that starts from static but unstable initial conditions. 
This could not be answered within the context in which the instabilities appeared, because 
that context was not a mathematically well defined model. Instead, Chandrasekhar's re- 
sults have been taken to mean that there are limitations to the range of physical parameters 
(mass and size) that are possible, given the assumed thermodynamic properties of the star. 
With this conclusion our results are in perfect agreement: There is a largest mass, beyond 
which the problem is not that the static solutions are unstable, but that they do not exist. 

The maximal mass of a neutron star can be obtained once the critical density is known. 
Using commonly accepted values we again recover the traditional limit. We explore some 
variations of the representation used for the piecewise almost constant function n[p] and 
discover an unexpected density profile. 
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2. Matter model and equations of motion 

We add the matter model action 

A matter = J ^X^^g^,^ - C 2 ) - V[p\) = : J d 4 X^9 C 

to the Einstein-Hilbert action and restrict the metric to the spherically symmetric form 

(ds) 2 = e u {dt) 2 - e x (dr) 2 - rW, 
Einstein's equations then reduce to 

G\ = -e" A (^ + y) + = 8nG(e-W a - c), 

G/ = -e" A + 1) + 1 = 8ttG( - e" A ^') 2 - £), (2.1 - 3) 

G\ = e - x ^ = -87rGe- A ^V, 
and are supplemented by the wave equations, from variation of the fields p and ip, 

Ve-^ 2 -e- A (^) 2 -l) = 

2 V yv> J dp' (2.4-5) 

a t ( e (-+A)/2 r 2^ ) _ (e (-A)/2 r 2 pV; /y = Q 

The Lagrangian density £, evaluated on-shell, is interpreted as the pressure, subsequently 
denoted p. This is not only because of its appearance in the expression for the energy- 
momentum tensor - in Eq.s (2.1-3) - but also because of the fact that C = p(dV/dp) — V, 
which is a familiar expression for pressure in terms of internal energy; see [FW], page 304. 
(Extending this to an off-shell identification of the pressure with the Lagrangian density 
would be a mistake.) 

The function A is often replaced by the function M defined by 

r 2M 
M:=-(l-e- A ), e" A = l-— ; (2.6) 



then Eq.s (2.1-2) can be written as follows, 

M / = 47rGr 2 (e-VV> 2 -£), 
re"V = 1 - e" A + 8nG r 2 (e~ x pip' 2 + C). 



(2.7) 



The two equations can be combined to yield 

{v + X)' = 8nG r e x p (e"> 2 + e" V 2 ). (2.8) 

The differences between this theory and the usual phenomenological one are mainly as 
follows: 
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• Eq.(2.4) looks unfamiliar, but taking the derivative with respect to r one recovers the 
usual force equation with only minor changes. In the static case, when ip = = and 
v, A are time independent, this becomes the hydrostatic condition, 

/ / 1 / —v 

V 19 = e 

and this has exactly the same form in both theories. But Eq.(2.4) is stronger than its 
derivative; it furnishes an additional constraint on the boundary, and it is this equation 
that provides a relation between the radius and the mass. In the traditional approach the 
boundary of a polytrope is often chosen to be at the point where the pressure becomes 
zero, that always exists if the polytropic index n = (7 — l) -1 is less than 5. There is 
no effective or meaningful matching of the field v to an external Schwarzschild metric, a 
fact that, in our opinion, makes the whole proceeding unsatisfactory. In our model this 
matching is expressed by 

u(R) + \(R) = 0, l- e "(K> = *£l. (2.9) 

-Tt 

The first equation determines R and the second gives the value of the mass. 

• Another significant difference is that the old approach has no intrinsic conserved current 
(a conserved baryonic current is often introduced by hand [C2]), while the new theory 
does, namely 

J» = V=gPirdvil>. (2.10) 

The existence of this conserved quantity does not by itself assure us that the mass is a 
constant of the motion. However, with the new equation of state introduced later, and 
the natural boundary conditions that it entails, it does indeed turn out that the mass is 
conserved. (Section 3.) 

• In the static case the pressure of the model corresponds exactly to the pressure as defined 
by Tolman's formula, while Tolman's density is replaced by p + (9V)'. We do not try to 
guess the precise physical interpretation of p and p but try instead to obtain results in 
terms of quantities that we feel sure are physical, such as the gravitational mass, uniquely 
defined by the asymptotic gravitational potential. 

In the static case, when tp = 1 and ifj' = 0, we find a curious, special solution, with / 
constant, namely 

e -^ = l + 2a7/, e" A = 1 + 8nGpr 2 . (2.11) 

As a global solution it is of no interest, but we shall find a solution for which / remains 
nearly constant over a finite interval. 

In the paper [F2] we examined static solutions with boundary conditions determined 
by matching the metric of the interior polytrope to an exterior Schwarzshild metric. As ex- 
plained above, the radius and the mass were determined by the two conditions in Eq.(2.9). 
Solutions were given for n = 1, 2, 3, 4, 6, 10, ... , but only the case n = 3 will be invoked 
here; Table 1 reproduces the data for this case. 



7 



Table 1. Mass/radius relation for n = 3 and exterior Schwartzschild. 
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There appears to be an upper limit to the mass of about 1.02, and a lower limit on the 
radius of about 1.96, in units where c = G = a = 1. The number in the third row is the 
position of the first zero of the pressure. Since the scale is not fixed the most significant 
data are the range of the dimensionless ratio 2mG/R and the correlation of this number 
with z/(0). This is shown in Fig.l, lower curve. 

Oscillations around the static solutions 

In this section we use units such that 8irG = 1 and fix the scale by setting a = 1. We 
linearize the equations at the static solution. The equation for G r r gives dp, 

re x 5p = bv' - [y' + -)SX, (2.11) 
r 

and so does "Newton's equation", Eq.(2.4) 

(e-^ 2 )5(e-^ 2 ^) = S^ = Up, 

dp p 

or 

5p = e - v p{5^ -Su/2). (2.12) 
Eliminating dp from these two gives 

rpe^idiP - 6u/2) = 6u' - {v' + -)6X. (2.13) 

r 

Equation (2.8) for G\ — G r r gives 8p, 

re x ~ v [dp + p5X + 2p(8ip - 6u/2)^ = 6u' + SX', 



or 



r5p = -2rp(5ip - Su/2) - rpSX + e v ~ x {5u' + 5X') (2.14) 
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Fig.l. The lower curve shows the essential information from Table 1. The abscissa is v(0) 
and the ordinate is the dimensionless ratio 2mG/R, where R is the boundary defined by 
Eq.(2.9). The upper curve shows the same information from Table 2, the ordinate is here 
2mG/R cr . 

Using (2.13) we eliminate the time derivative and obtain a constraint, 

r5p = e v ~ x (SX' - 5u' + (z/ - A' + 1)5X) (2.15) 



The original set of equations is equivalent to (2.12) and (2.15) that determine 5p and 
Sp, and the two equations 

rpe x - u (SiP - <W2) = Si/ - [y' + -)SX (2.16) 

r 

and 

SX = rpip' (2.17) 

for the functions A anf ip. The only difficulty is that there is no determination of du; to fix 
it we have to use the equation of state; the choice may affect stability. 

Equation of state 

We express V and p parametrically, in terms of the Emden function /: 

V = afp, p = f n , p = p^V, 

dp 

with 

n 2 -n 1 
n = n 1 + 1 + fK - 
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This function takes the constant value n\ (in most of our calculations n\ = 3) in the inner 
regions where / > 1, and the constant value ri2 (the values 6 and 15 were explored) in the 
atmosphere. The critical region is thus at the place where / takes the value unity; / and 
p are entered as multiples of their critical values. This choice, with a constant, gives the 
best approximation to an isentropic equation of state, as we shall see. 

In regions where n is constant we have p = (a/n)f n+1 . In the boundary region there 
will of course be deviations from this equation of state. We calculate 

dV 1 + n-Nlnf (n 1 -n)(n-n 2 ) „ 

dp n — JXImj n\—ni 

Because of the high value of K that was used (K = 50), the function N vanishes except 
in a narrow interval where ri\ < n < ri2- Of course, because K is large, iV is not small 
in this region. However, for the same reason, the function / varies little from its critical 
value in this interval. Therefore, taking this critical value to be unity, making / rs 1 in the 
interval, is highly beneficial. The two functions p-jj^V, *yV. are nearly indistinguishable 
over the entire interval < r < oo. The deviation is small, in a small region, and it will 
be ignored in the calculations. (Most results show a remarkable insensitivity to details of 
the equation of state near the critical point.) Thus we set p = ap 1 and 

pSp — -fpSp. (2-18) 

For a stationary solution, when SX = sSX for some number s that we hope will have to be 
pure imaginary, from (2. 11), (2. 15) and (2.18), 

Si/' = re x 5p+ iv' + -)SX, 
r 

rpip' = s SX, 

8X' = rKe x 5p + (A' - -)8X, K = 1 + e~ u php. 

r 

Elimination of 5p leads to 

6v' = re x ~ u p(s5ifj - 8v/2) + {y' + -)8X, 

r 

rpifj' = s SX, 

SX' = srpKe x - v SiP - {rpK /2)e x ~ u Su + (X' - -)SX. 

r 

or 

/ Su\ ( -{rp/2)e x ~ v s{rp)e x - v f ' + £ \ f Su\ 
(d/dr) \8ip = (s/rp) \ [ Sip . (2.19) 

\SX J \-(rpK/2)e x ~ v s(rp)Ke x ~ v X' - £ / \ SX J 

From this it is easy to see that the integration from r = can proceed, we start at r = 10 -10 
with Sip = SX = and Sv ^ 0. This will make Sv' and SX' of order r, SX of order r 2 and 
Sv — 5is(0) of order r 2 . 
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The solutions include a simple gauge transformation; when it is ignored the system can 
be reduced to a single, second order differential equation for the function L = re~ x SX: 

L = r {W + \)L + r^e-^^C—^)'. (2.20) 

This equation shows that the stationary solutions have frequencies determined by a self 
adjoint Sturm-Liouville operator, with a domain of functions L that satisfy a condition 
that fixes the value of V /L at the boundary and such that L/r 3 is regular at the origin. 

Boundary conditions, difficulties 

Round 1. Here we report the result of studying the stability of the static solutions in 
the case that the star has a fixed polytropic index in the interior and the metric is matched 
to an exterior Schwarzschild metric at the value of r determined by Eq.s (2.9). (These are 
the solutions listed in Table 1.) We adopt Eddington's boundary conditions at the center 
[E] , so that the solutions of the static equations of motion for fixed values of n are indexed 
by the value z/(0) of the function v at the center. Matching of the solution to the exterior 
metric - Eq.s (2.9) - determines both the radius and the gravitational mass. The result is 
a relation between mass and radius, for each value of n, reproduced for the case n = 3 in 
Table 1 and in Fig.l, lower curve. 

The stability of a star is to be determined by solving equations (2.19) or (2.20) with some 
boundary conditions. It is difficult, however, to understand what boundary conditions are 
appropriate. 

To deal with Eq.(2.20) from the point of view of Sturm-Liouville theory, one looks for 
solutions with harmonic time dependence, L(x,t) = e lujt L(x,0). One identifies the range 
of the parameter uj with the spectrum of an operator in a Hilbert space constructed from a 
space of functions of r. An acceptable set of boundary conditions must make this operator 
self adjoint. 

It has always been assumed that the origin is a regular point. It is difficult to find a 
real justification for this, since the center of the star is a region about which one has very 
little information. Nevertheless, we follow this precedent since it helps to give a precise 
mathematical sense to our model. Thus v has a well defined value at r = 0, A is of the 
order of r 2 and the function L is of order r 3 . In this case the possible additional boundary 
conditions amount to fixing the value of L'(r)/L(r) at the surface. If we fix the boundary 
and suppose that L = there, then we find a discrete set of oscillating solutions and, in 
some cases, decaying solutions, in agreement with the findings of Chandrasekhar [CI]. But 
the radius of the star is determined by the first of the conditions (2.9), and that means 
that the surface of the star is pulsating, as we have verified numerically. This seems not at 
all unnatural, the difficulty is that the asymptotic mass, as determined by the matching 
of the metrics, is also pulsating. We are very skeptical of these results. In the traditional 
treatment mass has no relation to a conserved quantity, and an uncertain relation to the 
asymptotic mass; this may be the reason why the problem has not been addressed in 
connection with the work of Chandrasekhar. 

Round 2. Fixing the mass would seem to be the more reasonable boundary condition, 
since it is determined asymptotically in a region where the density is zero. It does not 
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seem possible that oscillations of a finite star propagate to infinity through an infinite 
region of empty space. To understand this better we should remember that the external 
Schwarzschild metric, according to our interpretation [F2], is not the metric of empty 
space, but a singular limit of a family of metrics of spaces with non vanishing density. 
Consequently, in our next attempt we replaced the empty Schwarzschild exterior by a 
crude approximation for the atmosphere, another polytrope. The mass distribution now 
extends to infinity, though the density falls off extremely rapidly and the metric soon 
becomes indistinguishable from that of Schwarzschild. The asymptotic mass is a property 
of the exterior polytrope, but we were unable to match the two polytropes at the boundary 
in such a way that this mass would remain constant. Still we cannot claim to have excluded 
this possiblity completely. 

It has been traditional since the beginning, to admit a discontinuous behaviour of density 
and pressure at the surface of a star. This is reasonable if the star is cold, but perhaps 
less likely to be typical of a polytrope. Facing doubts of this kind, and the difficulties 
discussed in the preceding paragraph, we came to the realization that it may be better to 
give up the idea of a fixed boundary and introduce the equation of state described in the 
introduction, thus allowing all the fields to vary continuously throughout. 



3. Improved boundary conditions 

Round 3. We suppose that there is a region of critical density p cr , where the polytropic 
index changes more or less abruptly from a value n\ < 5 (in our calculations n\ < 3) that 
is appropriate for the bulk of the star, to a value > 5 (actually 6 or 15) that we hope 
may be appropriate for the atmosphere. Precisely, 

V[p\=ap cr p y , 7=1 + -, n= l + j K » W- 1 ) 

where K is a suitably large number (actually 50) and / = f / f cr , p = p/ Per are all dimen- 
sionless. The critical density now appears as a common factor of the energy momentum 
tensor, and on the right hand sides of Eq.(2.1-3), together with G. We shall drop the tildes 
on / and p, so that / and p are henceforth given as multiples of their critical values. Then 
Eq.s (2.1-5) remain valid if the factor G in (2.1-3) is replaced by Gp cr . Finally we choose 
our unit of length such that this factor is equal to unity, 

Gp cr = 1. 

The pressure in these units is (a/?i)/5 7 and the critical pressure is (a/n)p cr , n = {n\+ri2)/2. 

There is no longer any question of matching to an exterior Schwarzschild metric, instead 
we require that the metric approach the Schwarzschild form at large distances, to order 
1/r. The mass is determined by this asymptotic metric, 

2mG :— limrA(r) = — lim rv(r). 
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In the traditional approach the second condition is not effective, since only the derivative v' 
of the function v appears in Einstein's equation. The extra condition that comes from the 
action principle guarantees that the metric is asymptotically Schwarzschild (as r 3 p(r) — > 0) 
so that the two limits always coincide. The indices 711,712 were given the values 3,6 in our 
initial calculations, Table 2. A larger value of n,2 makes the metric approach more quickly 
to the Schwarzschild form, Table 3. The exponent K determines the abruptness of the 
change of the index from 3 to 6, our calculations were done with K = 50. 

Since the index is not constant, the value of a can no longer be reduced to unity by a 
change of scale. 

As before, we assume that all the fields are regular at the origin. A star is characterized 
by the parameters f cr and a. We choose a value of a and determine allowed value(s) of 
z/(0) by demanding that —ru(r) tend to a finite limit (twice the mass m) at infinity, and 
establish in this way a correlation between R and m. Results are given in Tables 2-7. 
The upper curve in Fig.l shows the dimensionless number 2mG/R cr versus f(0), for the 
case n = 3,6 (3 inside and 6 outside). The lower curve was obtained by matching the 
interior solution (n = 3) to an exterior Schwarzschild metric. The difference between the 
two curves is easily explained since the upper curve refers to the critical radius R cr while 
the lower curve is 2mG/R, where R is the boundary. 

All static solutions are stable! 

When there is no fixed boundary, and the star extends to infinity, the asymptotic 
behaviour becomes important; for the function L that embodies the oscillations around a 
static solution we find 



with k and b constant. Values of the exponent k determined by numerical calculations 
for n = 6 (k = 5/2) and n = 10 (k = 9/2) at large distances are such that the metric 
fluctuations 8v and 8\ fall of faster than 1/r. It follows that the mass (defined by the 
asymptotic metric) is unaffected by the oscillations. 

With the improved boundary conditions the real spectrum of frequencies is continuous 
and apparently the entire real line. No unstable solutions were found. This could have 
been anticipated by inspection of Eq.(2.20). The factor p in the first term on the right 
hand side makes this term fall off very fast at infinity, which suggests that this term cannot 
affect the spectrum of u 2 . The conclusion must be that the instabilities first discovered by 
Chandrasekhar [CI] are imposed on the theory by the choice of boundary conditions. 



L oc sin(r 3 / 2 6)/r fc 



Constants of the motion! 

The conservation law (2.5) can be integrated to yield 




In view of the boundary conditions at the origin, 




13 



The factor p on the right hand side suggests that there is no flux at infinity, but in fact 
the flux rptfj' is equal to — 5A/8-7T by Eq.(2.3). For a static configuration both sides of this 
equation are zero; for a first order deviation from a static configuration we have 



d_ 

dt 



I 



If the perturbed and unperturbed metrics both tend to Schwarzschild at infinity, then 



It is not a priori obvious that the asymptotic mass is a constant of the motion, but a 
result of our calculations is that r5X tends to zero at infinity so that in fact m = 0. The 
asymptotic mass is a constant of the motion and so is the quantity 



In our numerical study this number turns out to be bounded upwards by the asymptotic 
mass, the difference being greater in the case of strong gravitational fields, which suggests 
an interpretation of the difference as a binding energy. This is the interpretation suggested 
in [KW], in the context of the traditional treatment, with the difference that they replace 
p by p (Tolman's density). But the integral of p is not a constant of the motion! 



r8X — > 2m so that, finally, 




(3.2) 
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4. Applications 

Modelling the sun 

The results for n\ = 3, = 6 (Table 2) and for n\ = 3/2, -n-2 = 6 (Table 6) are shown in 
in Fig. 2. The graph of allowed values of (R cr , 2mG) has a lower branch with small internal 
pressure, a maximal value of R cr , and an upper branch with increasing central pressure 
and a maximum value of the mass. 
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Fig. 2. The mass-radius relation obtained with different boundary equations of state. On 
the left n = 3, 6, on the right n = 3/2, 6. 

The bulk of the sun is often modelled with a poly trope with index 3. The parameter 
values are 

(2MG) sun = 2.95 x 10 5 cm, R sun = 6.96 x 10 10 cm, 

and the ratio is .424 x 10~ 5 . We must not equate R sun with the value R cr of the radius 
at which the density takes the critical value; this may occur deep in the interior. 
If the critical density is p cr = k 2 g/cm 3 : then 

p cr G = .7414 x 10 - 28 k 2 /cm 2 

and the unit of length used in our calculations is therefore 



£ = 



1.16 



x 10 14 cm. 



k 



Here are 2 examples from Table 2. 



Example 1. Take 1/a = 10 5 , the Table gives 2mG = .32 x 10 -7 £ = 37/c -1 x 10 5 . To 
agree with the sun value we need k — 37/2.95 = 12.54. Thus p cr = 158 g/cm s , which may 
be a little high. The critical radius is 9.23 x 10~ 4 £ = .854 x 10 10 and the critical pressure 
is p cr = (a/3)p cr = .527 x 10" 3 or 4.74 x 10 17 dyn/cm 2 . 

Example 2. Take 1/a = 2 x 10 5 , the Table gives 2mG = .113 x 10" 7 £ = .131 k' 1 x 10 7 . 
Here we need k = 13.1/2.95 = 4.44. Thus p(0) = 1.257 x p cr = 24.8 g/cm s , which may be 
too low. The critical radius is 



R r 



6.52 x 10" 4 £ = 1.7 x 10 



10 
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This is about R sun /4. The density profile, shown in Fig. 3, shows that this is quite 
reasonable. The pressure is 

P(0) = ^Pcrf(0) 4 = .446 x 10" 4 , 

or .401 x 10 17 dyn/cm 2 . 




0.001 0.002 0.003 0.004 



Fig. 3. Various functions plotted against the distance from the center. The function with 
the step is n(r)/10. The other solid curves show f(r) and p(r). The dashed curves show 
analogous results with a polytropic index fixed at the value 3. What cannot be seen is that 
the latter pass to negative values at r = .0047, while the solid curves never do. 

Given the known values of the mass and the radius of the sun, the traditional polytrope 
model predicts the values p(0) = 76.39(7 /cm 3 and p(0) = 1.24 x 10 17 dyn/cm 3 for the 
density and pressure at the center. Our dynamical model gives slightly different values for 
the central density and pressure. Various refinements, such as larger values of K and/or 
ri2, may affect these predictions to a limited extent. 

Let us compare the effect of different boundary conditions for fixed values of u(0) and 
a. In the case of Example 2, in units of 10 10 cto: 

• The Emden function vanishes at 12.3. This would be the prediction for the radius of 
the standard approach with this choice of parameters (not the best choice). The density 
predicted by the model at this point is about 1/5000 of the critical value. 

• The actual value of the visual radius is 6.96. 

• Matching to an external Schwarzschild metric would fix the boundary at 3.9. The 
density at this point is about 1/8 of the critical value, dropping abruptly to zero. 

• The critical radius is at 1.7. 
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• The density profile is shown in Fig. 3, and the density predicted by the standard 
theory is shown for comparison. 

All three models can be made to give the correct radius and mass of the sun. They 
differ, but not very greatly, in the predictions for central density and pressure, and they 
differ considerably in the degree to which they seem to be physically reasonable. 

Fig. 4 shows profiles of the metric functions. 
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Fig. 4- The metric functions A (lower curve) and —v in the case that best fits the sun. The 
dashed line is the function —u for the Schwarzschild metric with the same mass. 

Strong gravitational fields 

A good indicator of the strength of the gravitational field is the dimensionless ratio 
2mG/R cr . The maximum values are, for the case that = 6: 

m = 3 5/2 2 3/2 1 1/10 
2mG 

—— = .33 .43 .50 .59 .67 .77 

As the value of is increased the falling off of the density outside the critical radius 
becomes much more rapid. Nevertheless, there is always a region where the gravitational 
field is less strong and where the density is not yet negligible. Thus it is difficult to 
imagine a situation where no radiation escapes. A much simpler assumption to explain 
the astrophysical "black holes" would be for the star, or at least the atmosphere, to be 
non-radiating on account of the temperature being very low. But in the polytropic models 
the temperature falls off much more slowly than the density. 
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The Chandrasekhar limit 

Under conditions that are believed to prevail in a white dwarf the value of the parameter 
a can be calculated. Normally the equation of state is expressed as p = Kp 1 , with pressure 
and density given in units of g/cm 3 . In that case one obtains 

K = .548 x 10" 6 , a = 3f cr K. 

Our unit of length is 

£ = x 10 14 cm. 

f6/I 
Jcr 

The numbers posted for 2mG, to be expressed in centimeters, must be multiplied by 
£, including the factor f C r 3 ^ 2 oc a -3 / 2 . To find the maximum value of the mass we have 
evaluated the product of 2mG in the table by this factor, with the result that it is essentially 
constant at 1.01 in the upper half of the table and eventually decreases to about .33 near 
the bottom of the table. The model is thus in agreement with the traditional treatment, 
in predicting a unique mass in the case of weak gravitational fields, and a maximum value 
of the mass that is reached roughly in the interval < a < .01. This maximum value is 

2mG£ = 1.16 x 1.01 x 10 14 (3K) 3/2 = M{2mG) sun . 

When the mass is close to this upper limit the ratio 

2mG/R cr < .334. 

The maximum value is 788 times the value of 2mG/R sun , so it would appear that a white 
dwarf of this mass may have a radius as little as 1/100 th that of the sun. 

The model makes these predictions, basically the same as the theory of Chandrasekhar, 
but without invoking instabilities. All the stars described by the model are stable. Field 
profiles are shown in Fig.s 3 and 4. 

Neutron stars 

According to Oppenheimer and Volkoff [OV], a model of a neutron star must have a 
very rigid equation of state. We tried n = 1/10 in the bulk and n = 6 outside and in this 
way we reach the higher value of .77 for the ratio 2mG/R cr . The model does not furnish 
the scale; that is, in our case, the value of p cr . It is usual to place the central density of 
a neutron star between 10 14 and 10 15 g/cm 2 . Taking p cr = 10 14 g/cm 2 gives us the scale 
£ = 1.16 x 10 7 . Two sample solutions are 

a = 1/2, 2mG = 9 x 10 5 cm, R cr = 17 km, 

which is 3 solar masses and 2.4 x 10~ 5 solar radii and 

a = 4/3, 2mG = 10 6 , R cr = 13 km. 

If the density is 9 times greater the numbers will be 3 times smaller. All these numbers 
are consistent with current estimates. 
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It was found that the formula that was used for the variable poly tropic index fails, in 
this case of a low internal value, to give a constant value in the interior; the value is close 
to n\ near the center, but begins to increase at about half the critical radius. The natural 
remedy would be to increase the value of the exponent K, but in that case the power of 
Mathematica to handle very large numbers is overtaxed. An alternative formula, namely 

m+n 2 ni-n 2 / - 1 
n\j\ = ; = (4.1) 

approaches a step function very well if e is very small, and for 3 > ri\ > 1 this formula 
reproduces the same relation between mass and radius as the one used to construct the 
tables, to a very good accuracy. However, the density profile is affected in an interesting 
way as e is decreased. In an intermediary range of the radial variable the solution slips 
into the solution with constant density described in Section 2. In this interval p ~ p cr and 
the polytropic index is determined by the metric, more precisely by z/, through equation 
(2.11). This phenomenon is specific to the model, to the use of a variational principle. 
Consequently, we succeed in fixing n = n\ in the region interior to this platform only. 

As shown in Fig. 5, the effect is not very important in the case ri\ = 3, n 2 = 6, but for 
smaller values of ni it is decisive. Fig. 6 shows the density profile for the case ri\ = 2, n 2 = 6 
and e = 10" 6 . 

Thus it seems that attempts to fix n at a very low value leads instead to a stratification 
of the star, with a middle region in which the density is constant. This has an uncanny 
similarity with the accepted picture of neutron stars. If the value of ni is decreased below 
.3 or .2 this middle region reaches the center, and the theoretically uncertain core shrinks 
to nothing. We do not have the temerity to suggest that this picture of a neutron star 
corresponds to reality, but we find it fascinating. 



4. Discussion of the results 

1. The success of the polytropic equation of state in accounting for a wide range of 
astronomical objects is well known and almost miraculous. One may feel, nevertheless, 
that there is some question about the best choice of boundary conditions. There is also 
something a little unsatisfactory about the definition of mass and mass density. We have 
advocated the use of a variational principle, and we have found that there is a simple and 
natural matter lagrangian that allows us to reproduce all of the results of the phenomeno- 
logical approach. It also provides a conserved quantity, something that is often added as 
an additional ingredient to Tolmans theory. The new equation of state, that interpolates 
between an interior polytrope and an exterior polytrope that extends to infinity, is justi- 
fied on physical grounds; see the Introduction. It wipes out all the uncertainties that are 
presented by the more traditional approach using a fixed boundary [E][C1], and it allows 
to establish the existence of a constant of the motion related to the conserved current. 
The idea of a double polytrope is not new, for Chandrasekhar has proposed a very similar 
equation of state [C2],[OV]. With this equation of state all the examined double polytropes 
are stable, which is an extremely satisfactory resolution of the paradox presented by Chan- 
drasekhar's limit. It is worth emphasizing that the only practical test of Chandrasekhars 
prediction is the observation of an upper limit on the mass; this is a prediction of our 
model as well. 
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Fig. 5. Contrast between two ways of smoothing the transition at the critical density. On 
the left the formula (3.1), on the right ( 4-1)- ^ n both cases, the parameters used are those 
that best fit the sun. 

2. In this paper it is taken for granted that the metric approaches the Schwarzschild 
form at infinity; in particular that the function rv(r) has a limit at infinity. What is called 
"the mass" , and denoted m, is defined by 



The traditional approach defines "the mass" as the integral (in traditional notation) 



its value is the same as the asymptotic mass. Of all the textbooks that we have consulted 
only one ([KW], pages 12-13) expresses any discomfiture with this formula; they say that it 
is "treacherous" , because the measure is wrong, given that p is a scalar field. The difficulty 
in the traditional theory arises because it does not have a conserved current, ultimately 
because it is not formulated as an action principle. Our model is, and we have shown that, 
with the boundary conditions that we have adopted, there is a conserved quantity, 



The value is always less than m and the difference can perhaps be interpreted as a binding 
energy, see [KW] page 387. To this we must say that the only meaningful concept of 
energy in General Relativity is the ADM energy; therefore we make no suggestion as 
to the proper name for the above integral. The important point is that both it and the 
asymptotic mass are constants of the motion. It is our opinion that this solves a difficulty in 
the applications of General Relativity to mass distributions in general and to the structure 
of stars in particular. 



2mG = — lim ru(r). 



r—>oo 
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Fig. 6. Here the formula (4-1) was used with n = 2, 6. The platform of constant density 
reaches almost to the center. The lower curve is the function n(r)/10. 

2. The paradox that is presented by the Chandrasekhar limit is that his theory fails 
to explain what happens to an unstable configuration. Our model reproduces the limits, 
but the unstable configurations of Chandrasekhar's theory do not have counterparts in the 
model. The problem of explaining or predicting what happens to them simply does not 
arise. 

3. Observation of most real stars depends on the fact that they are luminous. The role 
of radiation in determining the equation of state has been overlooked in this paper and this 
is an obstacle to further phenomenological applications. The question of how radiation is 
to be included in the model is interesting. Of course, matter has to interact with radiation, 
but since the interior of a star is neutral on a large scale this may not be the first item to 
take up. Radiation can be introduced by adding the Maxwell action to the gravitational 
and matter actions, but this is probably not the best way. Here we suggest, but offer no 
justification for it, to include a contribution to the action of the form 



The modification of the Maxwell action by the factor a and the internal energy W is 
necessary for the same reason that we need the factor p and the internal energy V in the 
matter action. It is probable that the gravitational action should be similarly modified, to 
reflect the existence of a background of soft gravitons. Compare [R] and [AMM] . 
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where a is a scalar field, the density of radiation, and 




21 



References 

[AMM] Antoniadis, L, Mazur, P. and Mottola, E., 

Cosmological Dark Energy: Prospects for a Dynamical Theory, |gr-qc/0612068 
[CI] Chandrasekhar, S., The maximum mass of ideal white dwarfs, 

Astrophys. J. 74, 81- (1931). 
[C2] Chandrasekhar, S., Stellar Configurations with Degenerate Cores, 

Monthly Notices R.A. 95, 226-260 (1935). 
[Ed] Eddington, A.S., The Internal Constitution of Stars, Dover, N.Y. 1959. 
[Em] Emden, R., Gaskugeln, Teubner 1907. 
[Fl] Fr0nsdal, C, Growth of a Black Hole, 

J.Geom.Phys. 57, 167-176 (2006). [kr-qc/0508048| vl 
[F2] Fronsdal, C, Ideal Stars and General Relativity, |gr-qc/0606027 



[H] Hartle, J.B., Bounds on mass and moment of inertia of non-rotating 

stars, Physics Reports, 46, 201-247 (1978). 
[KH] Kippenhahn, R. and Weigert, A, "Stellar Structure and Evolution" , 

Springer- Ver lag 1990. 
[ MM] Mazur, P. and Mottola, E., Gravitational Vacuum Condensate Stars, 

^q"c/ 0407075 



[OV] Oppenheimer, J.R. and Volkoff, G.M., On Massive Neutron Stars, 

Phys. Rev. 55, 347-381, (1939). 
[R] Rees, M.F., Effects of very long wavelength primordial gravitational radiation, 

Mon.Not.astr.Soc. 154 187-195 (1971). 



22 



v(0) (l/o) 2mG [R cr ] 



.3469775 10 5 


(825 000) 


1.35 10~ 


9 


[3.22 10~ 


-4 


_ 5 

.4089376 10 


(700 000) 


1.72 10~ 


9 


[3.47 10~ 


-4 


.477093666 10~ 5 


(600 000) 


2.175 10 


-9 


[3.77 10~ 


-4 


.5725120 10 


(500 000) 


2.83 10~ 


9 


[4.13 10~ 


-4 


.715640106 10~ 5 


(400 000) 


4.00 10~ 


9 


[4.62 10~ 


-4 


.95419005 10 -5 


(300 000) 


6.155 10 


-9 


[5.3410 _ 


4 ] 


.143128161 10 _ 4 


(200 000) 


1.13 10~ 


8 


[6.52 10~ 


-4 


.19083741 10~ 4 


(150 000) 


1.74 10~ 


8 


[7.58 10~ 


-4 


.28625595 10~ 4 


(100 000) 


3.2 10 


i 


[9.23 10~ 


-4 


.572512 10 _ 4 


(50 000) 


.905 10~ 


7 


[.000585] 




.00057250563 


(5 000) 


2.9 10 


[.00415] 




.0023853416 


(1 200) 


.00002425 


[.00844] 




.0140868 


(200) 


.000349 




[.020] 




.02861153 


(100) 


.00097 




[.029] 




.0476807 


(60) 


.00202 




[.038] 




.0572184 


(50) 


.00261 




[.042] 




.071534195 


(40) 


.00357 




[.0465] 




.0954368 


(30) 


.00529 




[.0535] 




.1435468 


(20) 


.00905 




[.0655] 




.1692845 


(17) 


.01110 




[.0687] 




.18602925 


(15.5) 


.01244 




[.0736] 




.1923886 


(15) 


.01295 




[.0748] 




.2035445 


(14.2) 


.013825 




[.0766] 




.2230195 


(13) 


.01534 




[.0796] 




.242466 


(12) 


.01680 




[.0825] 




.2658189 


(11) 


.01853 




[.0855] 




.2944944 


(10) 


.02053 




[.0884] 




.330789 


(9) 


.02286 




[.0920] 




.3788565 


(8) 


.02550 




[.0950] 




.447908 


(7) 


.02846 




[.0980] 




.575115 


(6) 


.03060 




[.0962] 




.67760 


(5.75857) 


.0290 




[.0910] 




.7529185 


(5.9) 


.02596 




[.0840] 




.790723 


(6.1) 


.02390 




[.0793] 




.874 


(7) 


.01833 




[.0680] 




.92404 


(8) 


.01466 




[.0605] 




.98369 


(10) 


.01039 




[.0512] 




1.00924 


(11.245) 


.00874 




[.0474] 




1.013842 


(11.5) 


.00846 




[.0468] 




1.022385 


(12) 


.00795 




[.0455] 




1.051813 


(14) 


.0064 




[.042 ] 




1.119438 


(20) 


.00395 




[.0353] 




1.17798 


(25) 


.00305 




[.0332] 




1.19458 


(26) 


.00295 




[.0331] 




1.2183451 


(27) 


.0028 




[.0332] 





/(0) 2mG/a a ' 2 

1.07918 1.011 

1.07918 1.007 

1.07918 1.011 

1.07918 1.006 

1.07918 1.012 

1.07919 l.on Table 2 

1.07919 1.011 

1.07919 1.011 

1.07919 1.012 

1.07921 1.012 

1.07942 1.025 

1.08017 1.008 

1.08526 1.0115 

1.09174 .97 

1.10101 .939 

1.10591 .923 

1.11359 .903 

1.17942 .869 

1.1579 .809 

1.17601 .778 

1.18846 .759 

1.19334 .752 

1.20209 .740 

1.218 .719 

1.23475 .698 

1.25606 .676 

1.28418 .649 

1.3232 .617 

1.3818 .577 

1.4832 .527 

1.749 .450 

2.09284 .401 

2.48505 .372 

2.7564 .360 

3.6683 .339 

4.5583 .332 

6.2787 .329 

7.3522 .330 

7.5735 .330 

8.0092 .330 

9.780 .335 

14.473 .353 

21.073 .381 

22.446 .391 

24.114 .393 



23 



Table 3, n = 3, 15 



"(0) 

-6 

5.821298 10 
3.146561 10 -4 
.0314001 


(l/o) 
(540 000) 
(10 000) 
(100) 


2mG 

— 9 

2.32 10 
-7 

9.2 10 

-4 

8.82 10 


[flcr] 

[5.5 10 

[.004] 

[.04] 


.06271548 


(50) 


.002387 


[.053] 


.15696 


(20) 


.00828 


[.08] 


.3221575 


(10) 


.01878 


[■1] 


.4165655 


(8) 


0.2326 


[.103] 


.4980215 


(7) 


.02572 


[.104] 


.64313 


(6.2) 


.0262 


[■1] 


.787609 


(6.5) 


.02135 


[.08] 


.965793 


(10) 


.01026 


[.055] 


1.1022975 


(15) 


.00386 


[.038] 


1.150818 


(25) 


.00292 


[.035] 



fcr 

1.644 

1.64887 

1.712 

1.7876 

2.0708 

2.896 

3.7255 

4.864 

9.236 

24.91 

227.0 

3431.4 

8312.7 



Table 4, n = 2.5, 6 



i/(0) 


(1/a) 


2mG 




fcr 


0.5582799 


(50.5) 


.0026 


[.0401] 


1.0578 


.0560465 


(50.3) 


.00261 


[.0401] 


1.0578 


.056381205 


(50) 


.002642 


[.0403] 


1.0579 




(35) 


.00435 


[.0481] 


1.0637 


.1399534 


(20) 


.00921 


[.064] 


1.0813 


.231989 


(12) 


.01735 


[.0815] 


1.1205 


.309115 


(9) 


.02400 


[.0924] 


1.1645 


.3995988 


(7) 


.0310 


[.1012] 


1.2281 


.4708388 


(6) 


.0357 


[.1068] 


1.2886 


.580451 


(5) 


.0411 


[.111] 


1.4051 


.835304 


(4) 


.0442 


[.108] 


1.865 




(3.95) 


.04318 


[.1053] 


8.152 


.9446 


(3.95) 


.0414 


[.1016] 


2.2173 


.99335 


(4) 


.03946 


[.092] 


2.429 


1.213762 


(5) 


.0267 


[.0775] 


4.2252 


1.309025 


(6) 


.02094 


[.069] 


5.7912 


1.380006 


(7) 


.01735 


[.064] 


7.437 


1.412063 


(7.5) 


.01601 


[.0621] 


8.315 


1.449928 


(8.1) 


.01472 


[.0605] 


9.439 


1.510485 


(9) 


.01424 


[.0595] 


11.34 


1.55061 


(9.5) 


.012625 


[.059] 


12.60 


1.60489 


(10) 


.01218 


[.0586] 


14.20 


1.62 


(10.1) 


.01213 


[.0588] 


14.62 
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1/(0) 


(l/o) 
( ) 


2mG 
.02127 


[flcr] 

[.085] 


fcr 

1.06765 








.272783 


(10) 


.0217 


[.0866] 


1.06864 








.286664 


(9.5) 


.023 


[.0886] 


1.07172 








.3020432 


(9) 


.02441 


[.0908] 


1.0754 








.3384081 


(8) 


.0277 


[.0955] 


1.08537 


Table 5, 


n = 


2,6 


.3601397 


(7.5) 


.02967 


[.0982] 


1.09232 








.3849225 


(7) 


.03184 


[.1008] 


1.10127 








.4134867 


(6.5) 


.03425 


[.1037] 


1.11298 








.4468366 


(6) 


.037 


[.107] 


1.12849 








.5343378 


(5) 


.04356 


[.1135] 


1.17742 








.671471 


(4) 


.0518 


[.120] 


1.2762 








.9584545 


(3) 


.05972 


[.122] 


1.6077 








1.15677 


(2.78) 


.0579 


[.1154] 


2.0198 








1.24684 


(2.78) 


.0552 


[.111] 


2.2975 








1.434277 


(3) 


.0471 


[.0995] 


3.1966 








1.7143226 


(4) 


.0333 


[.0833] 


6.071 








1.91458 


(5) 


.02714 


[.0781] 


9.6402 








2.03779 


(5.45) 


.02569 


[.078] 


12.124 








"(0) 


(1/a) 
(50) 


2mG 
.00266 


[ficr] 
[.0367] 


fcr 

1.02525 








.0925715 


(30) 


.00542 


[.0481] 


1.0279 








.11082561 


(25) 


.00696 


[.053] 


1.02922 








.18293098 


(15) 


.01347 


[.069] 


1.03453 








.2272555 


(12) 


.0177 


[.0772] 


1.03793 








.2472169 


(11) 


.0196 


[.08055] 


1.03951 








.271095 


(10) 


.0218 


[.0843] 


1.04145 








.2998917 


(9) 


.0246 


[.0885] 


1.0439 


Table 6, 


n = 


1.5,6 


.335651 


(8) 


.0279 


[.0933] 


1.0471 








.3811072 


(7) 


.0321 


[.099] 


1.05151 








.4088025 


(6.5) 


.0346 


[.102] 


1.05444 








.4408607 


(6) 


.0374 


[.105] 


1.05813 








.478416 


(5.5) 


.0406 


[.1085] 


1.06296 








.5230489 


(5) 


.04425 


[.1124] 


1.06966 








.5770453 


(4.5) 


.0484 


[.1164] 


1.07976 








.643911 


(4) 


.0532 


[.1205] 


1.09704 








.729582 


(3.5) 


.0587 


[.1245] 


1.1308 








.8457205 


(3) 


.0649 


[.1285] 


1.19687 








1.02041 


(2.5) 


.0712 


[.1308] 


1.33075 








1.40847 


(2) 


.0735 


[.1251] 


1.85382 








1.715 


(2) 


.0658 


[.1132] 


2.734 








2.112565 


(2.5) 


.05086 


[.0977] 


5.452 








2.358166 


(3) 


.0441 


[.0928] 


8.6144 








2.467983 


(3.2) 


.0424 


[.0923] 


10.367 
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Table 7, n = 1 



"(0) 


(l/o) 


2mG 


[Rcr] 


fcr 


.02791529 


(100) 


.000991 


[.0237] 


1.01622 


.11050725 


(25) 


.00698 


[.0501] 


1.02071 




(20) 


.00938 


[.0566] 


1.02212 


.1823652 


(15) 


.01353 


[.066] 


1.02441 


.2700904 


(10) 


.02197 


[.0814] 


1.02883 


.2988242 


(9) 


.0247 


[.0852] 


1.03027 


.4388956 


(6) 


.03766 


[.1026] 


1.03747 


.5202716 


(5) 


.04457 


[.1101] 


1.04193 


.639032 


(4) 


.0537 


[.1187] 


1.04909 


.721724 


(3.5) 


.05936 


[.1232] 


1.05483 


.8296405 


(3) 


.0658 


[.1278] 


1.0639 


.9772556 


(2.5) 


.0733 


[.132] 


1.08228 


1.1965135 


(2) 


.08135 


[.1346] 


1.1563 


1.630761 


(1.5) 


.0869 


[.1314] 


1.54041 


1.842675 


(1.4) 


.0854 


[.1268] 


1.85969 


2.250955 


(1.4) 


.0773 


[.1157] 


2.79388 


2.47326 


(1.5) 


.07144 


[.1098] 


4.07289 


2.854157 


(1.8) 


.0625 


[.1032] 


7.36191 


2.966246001 


(1.9) 


.0607 


[.1022] 


8.74897 


3.0853946 


(2) 


.0592 


[.102] 


10.438 
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